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ABSTRACT 

We investigate the plausibility of detecting X-ray emission from a stellar jet that impacts 
against a dense molecular cloud. This scenario may be usual for classical T Tauri stars with 
jets in dense star-forming complexes. We first model the impact of a jet against a dense cloud 
by 2D axisymmetric hydrodynamic simulations, exploring different configurations of the ambient 
environment. Then, we compare our results with XMM-Newton observations of the Herbig-Haro 
object HH 248, where extended X-ray emission aligned with the optical knots is detected at the 
edge of the nearby IC 434 cloud. Our simulations show that a jet can produce plasma with 
temperatures up to 10^ K, consistent with production of X-ray emission, after impacting a dense 
cloud. We find that jets denser than the ambient medium but less dense than the cloud produce 
detectable X-ray emission only at the impact onto the cloud. From the exploration of the model 
parameter space, we constrain the physical conditions (jet density and velocity, cloud density) 
that reproduce well the intrinsic luminosity and emission measure of the X-ray source possibly 
associated with HH 248. Thus, we suggest that the extended X-ray source close to HH 248 
corresponds to the jet impacting on a dense cloud. 

Subject headings: hydrodynamics - Herbig-Haro objects - ISM: jets outflows - ISM: individual objects: 
HH 248 - X-rays: ISM 


1. Introduction 

Classical T Tauri stars (CTTSs) are character¬ 
ized by being surrounded by a disk of gas as a 
result of the conservation of angular momentum 
during the star-formation process. The central 
star accretes material from the disk through mag- 
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netic funnels (|Koenigll 119911) . In addition, dense 
gas from the inner region of the disk is collimated 
into a jet as explained in the context of the widely 
acce pted theory of mag neto-centrifugal launch¬ 
ing ( Ferreira et al.l l2QQ6f) . Recently, magneto¬ 
hydrodynamic (MHD) simulations have been per¬ 
formed to evaluate the importance of the accretion 
disk magn etic field in the launching process of a 
stellar iet (iMatsakos et al. 2008 . 20091: Stnte et al 


2014t IStaff et ^l2014l) . The role of the magnetic 
field for the collimation of stellar jets in laboratory 
plasma experiments ha s been recently explored by 
Albertazzi et al. ( 20141) through detailed numeri¬ 
cal simulations and comparison with observations. 


When escaping from the stellar system, the 
jet moves through circumstellar and interstellar 
material, interacting with the ambient medium 
and producing shocks. Usually, stellar jets are 
revealed in the optical by the presence of a 
chain of knots, the so-called Herbig-Haro (HH) 
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objects. These knots are known to be associ¬ 
ated with the jet’s shock front and post-shock 
regions. Supersonic shock fronts and post¬ 
shock regions along the jet are detected in a 
wide wavelength range, from radio to optical 


band (|Reipurth fc Ballv 2001 . and references 
therein) and may be detected also in the ul¬ 


traviolet 




Gomez de Castro fc von Rekowski 


I2OIII: ICoffev et~aL I2OI2I : Schneider et al. 2013I ). 


The knotty structure of HH objects within the jet 
axis is interpreted as the consequence of the puls¬ 
ing nature of the ejectio n of material by the star 
(see lBonito et al.lboiOal ). A possible explanation 
for this pulsing nature is the variable nature of 


the stellar wind ( Matsakos et al.l 20091 ). whether 


through stellar magnetic cycles or by variations in 
the accretion rates. 

Pravdo et al.l ( 200l[ ) observed high energy emis¬ 


sion from high speed HH jets. In particular, hy¬ 
drodynamic (HD) models predict X-ray emission 
from mechanical heating due to shocks produced 
by the int eraction between th e jet and the ambient 
medium ( Bonito et al.ll20o3) . This is particularly 
true whe n the jet is less de nse than the ambient 
medium (|Bonito et al.ll200^ . Instead, X-ray emis¬ 
sion from heavy jets, i.e. jets that are denser than 
the ambient, is weaker. The analysis of X-ray data 
from stellar jets can yield constraints to initial con¬ 
ditions of the process such as the jet velocity and 
density. However, the low statistics in the few stel¬ 
lar jets detected in X-rays s o far makes it diffi cult 
to achieve robust results (see Giidel et al.ll2QQ5 . for 
an example). 

X-ray emi ssion from pro t ostella r jets was firstly 
reported by IPravdo et al.l ( 2QQl[ ). who detected 
an X-ray source coincident with HH 2 that was 
not associated with any other galactic o r extra- 


galactic sour c e (see IVelnsamv et al.ll2QI4j ). Later, 


Favata et al.l (j2QQ2f ) detected X-ray emission from 
HH 154, a jet associated with the protostar 
LI551 IRS5 in the Taurus Molecular Glou d (see 
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also iBallv et al.l I2QQ3I: iBonito et al 


S). 


20111 


Since then, the number of 


increased: HH 80/81 (IPravdo et al. 

20041. TKH 8 

(iTsiiiimoto et, a1.ll2nn4. HH 540 ( 

Kastner et al. 

200511. HH 210 ( 

iGrosso et al. 20061). HH 216 

(|Linskv et al.ll2QQ71 

1), HH 168 (|Pravdo et al.ll2009[). 


Several jets from other young st ellar objects 


^ reported lik ely detection from a jet from the 
errs DG rau fsee lGiidel etHll2nn8l:ISchneider et, a1 


20131 for recent results). Skinner et al.l ( 2011 1 de- 


tected evidence for extended X-ray structure in 
RY Tan. Finally, X-ray emission from a jet aris¬ 
ing from the FU Ori t y pe st ar Z GMa was de¬ 
tected by IStelzer et al. (|2QQ9[ ). A list of X-ray 


(YSO) were detected in X-rays, too. iGhdel et al 


de tections and propert ies of stellar jets is given 
in iBonito et al.l (2010b). The analysis of the X- 
ray spe ctrum of stellar jets divides them into two 
classes ( Bonito et al.ll20I0bl ): (I) X-ray sources de¬ 
tected at the base of the jet (within 2000 AU from 
the stellar source), with temperature > 2 x 10^ K 
and (2) X-ray sources detected at large distances 
from the jet source (several thousands of AU), 
characterized by high luminosity and low temper¬ 
ature. The interpretation of this behavior is still 
not clear, but the difference between both classes 
seems to be related to the location of the X-ray 
emission. 

In this work, we study a different scenario for 
X-ray emission observed from a jet: a stellar jet 
that moves through the interstellar medium (ISM) 
and impacts a molecular cloud (wall) with a den¬ 
sity several times higher than that of the ambient. 
In this scenario, the jet is heavier than the ISM 
but lighter than the wall. We perform simulations 
with different initial conditions, widely exploring 
the parameter space. To test our results, we com¬ 
pared them with the jet HH 248, detected in X- 
rays with the XMM-New ton mission (Src. 12 in 
Lopez-Garefa et al.ll2QI3[ ). 

The organization of this article is as follows. In 
Section [21 we present the numerical setup for our 
simulations. Results for the different initial condi¬ 
tions explored by us are discussed in Section |3l We 
then study the case of HH 248. A brief description 
of the system is carried out in Sections jH Specific 
simulations for HH 248 are performed in Section |5l 
Final discussion and conclusions are presented in 
Section |6l 

2. Hydrodynamic model and numerical 
setup 

Our model describes the impact of a proto¬ 
stellar jet against a dense molecular cloud. We 
adopted the hydro dynamic model discussed by 
Bonito et al. I (|2QQ7l hereafter BOP07), extended 
to include an ambient environment characterized 
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Fig. 1.— Schematic of the simulation domain 
setup. Boundary condition at the left border is 
reflective/axisymmetric. Outflow conditions are 
permitted in the remaining boundaries. 


Dense cloud 


ISM/ambient 


Jet 


by the presence of a dense cloud. The jet im¬ 
pact is modeled by numerically solving the time- 
dependent hydrodynamic equations of mass, mo¬ 
mentum, and energy conservation in a 2D cylindri¬ 
cal coordinate system (r, ^), assuming axisymme- 
try (see BOP07 for details). The model includes 
the radiative losses from optically thin plasma, 
and the thermal conduction, including the effects 
of heat flux saturation. The model is implem ented 
using the FLASH code ( Frvxell et al with 

the Piecewise Parabolic Method (PPM) solver 
that is better adapted for compress ible flows with 
shocks ( Colella fc Woodwardlll984l ). 

The jet is assumed to propagate along the sym¬ 
metry axis, z, through an unperturbed ambient 
medium with a dense molecular cloud. The cloud 
is described by a higher density region (wall) in 
pressure equilibrium with its surrounding and is 
situated at ^ 9000 AU from the base of the do¬ 
main (see Figure [T]). We chose this distance as 


we want to compare the results from the simula¬ 
tions with the observations of HH 248 (see Sec¬ 
tion [5] for more details). The computational do¬ 
main extends ^ 1400 AU in the r direction and 
^ 14000 AU in the 2 : direction. At the coarsest 
resolution, the adaptive mesh al gorithm used in 
the F LASH code (PARAMESH: iMacNeice et al. 


2000[ ) uniformly covers the 2D computational do¬ 


main with an initial mesh of 1 x 10 blocks, each 
with 8^ cells. We allow for five levels of refine¬ 
ment, with resolution increasing twice at each re- 
finement level . The refinement criterion adopted 
(|LQhnerll2008[ ) follows the changes in density and 
temperature. This grid configuration yields an ef¬ 
fective resolution of ^ 7 AU at the finest level, 
corresponding to an equivalent uniform mesh of 
128 X 1280 grid points. 

We explore the parameter space defined by 
the density ratio between the ISM and the jet 
(^ambient/'^jet): the density ratio between the wall 
and jet (nwaii/^jet)? and the jet’s velocity (ujet)- 
The aim is to investigate which physical condi¬ 
tions result in the production of detectable X-rays 
from the stellar jet afte r impacting the dense cloud 
(preliminary results in lOrellana et al. 2012 b The 
grid of explored parameters is chosen in order to 
explore in detail the case of a jet over-dense in 
the ISM that becomes under-dense in the cloud. 


Acco rding to previous results (e.g. iBonito et al 


I 2 OO 7 I ). this is the most promising scenario for X- 
ray emission. For completeness, we also explore 
the remaining two possible scenarios in less de¬ 
tail: the jet being heavier than both the ambi¬ 
ent and the molecular cloud and the jet being 
lighter than both the ambient and the molecu¬ 
lar cloud. The initial temperature of the jet is 
fixed to the value Tjet = 10^ K in all the simu¬ 
lations. Figure [2] shows the parameter space for 
our simulation. The shaded area represents the 
region of the parameter space in which the am¬ 
bient would be denser than the molecular cloud. 
This scenario is not explored in this work. In 
BOP07, it is shown that detectable X-ray emis¬ 
sion originates from the interaction of light jets 
with the ambient medium. Here we explore the 
scenario of a jet heavier than the ambient but 
lighter than the molecular cloud. This scenario 
corresponds to the top-left quadrant in the figure. 
Table [T] summarizes the parameters of our simu¬ 
lations. The size of the circles in Figure [2] cor- 
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Fig. 2.— Parameter space of the simulations per¬ 
formed in this work. Small circles correspond to 
jet’s velocities ujet < 500 kms“^, medium-size cir¬ 
cles are for 500 < ujet < 1200 kms“^, and large 
circles are for simulations with 1200 < ujet < 
2000 kms“^ and ujet > 2000 kms“^, respectively. 
The shaded area is the scenario with a molecular 
cloud less dense than the ambient, which is not 
investigated in this article. The continuous line 
corresponds to nwall/^ambient = 100 . 

responds to four different velocity ranges for the 
jet: Ujet < 500 kms“^, 500 < ujet < 1200 kms“^, 
1200 < Ujet < 2000 km s“^ andujet > 2000 km s“^, 
respectively from smaller to larger circles. In most 
cases, we simulated jets with different initial ve¬ 
locities for the same values of riambient/^jet and 

^wall/^jet • 

3. Results from simulations 

The evolution of the density and temperature 
of a jet with distinct ambient-to-jet density ra¬ 
tios (nambient/^jet) when traveling through a ho¬ 
mogeneous ambient was investigated in detail in 
BOP07. In such work, it was shown that the case 
of a light jet (nambient/^jet > 1) with high Mach 
number (Mjet) is the most favorable for detection 
of X-rays. A hot (T > 10^ K) and dense blob is 
formed behind the front shock in both light and 
heavy (nambient/'^jet < 1) j^ts whcn the jet’s ve- 


Table 1: Summary of the simulations performed in 
this work. 


^ambient /^jet 

^wall/^jet 

Mach 

'^jet ^ 

logTmax^ 




(kms-^) 

( K ) 

0.1 

0.5 

20 

936 

6.70 

” 

0.5 

30 

1404 

6.84 

” 

2.0 

20 

936 

6.62 

” 

2.0 

30 

1404 

6.80 

” 

10.0 

20 

936 

6.57 

” 

10.0 

30 

1404 

6.84 

” 

20.0 

10 

448 

5.70 

” 

20.0 

20 

936 

6.54 

” 

20.0 

30 

1404 

6.83 

” 

40.0 

20 

936 

6.64 

” 

300.0 

150 

2000 

7.33 

0.2 

20.0 

20 

662 

6.04 

0.5 

20.0 

40 

800 

6.41 

1.0 

10.0 

20 

296 

4.70 

” 

10.0 

200 

2960 

7.18 

” 

40.0 

100 

1480 

6.95 

” 

100.0 

100 

1480 

7.08 

10.0 

20.0 

100 

448 

5.10 

” 

100.0 

100 

448 

5.60 

” 

100.0 

300 

1404 

6.91 

100.0 

1000.0 

100 

148 

4.40 

- 

1000.0 

1000 

1480 

7.20 


^Derived from the Mach number. 

^ Measured inside the dense cloud, at 2 ; ~ 13000 AU (see 
Figure O. 

locity is high enough, but in heavy jets this blob 
is fainter in X-rays and considerably less extended 
(see Figure 8 in BOP07). In fact, the shock ve¬ 
locity required to produce this faint X-ray emis¬ 
sion in heavy jets is too high compared to obser¬ 
vations (ush > 2000 kms“^). At lower velocities, 
the X-ray luminosity of the heavy jet is low (see 
BOP07). Consequently, detection of X-ray emis¬ 
sion from heavy jets is not expected. Instead, the 
jet’s luminosity increases rapidly when it becomes 
lighter than the ambient (see Section O Figure [7]). 

To illustrate the previous discussion, we com¬ 
pare the evolution of the mass density and temper¬ 
ature of a heavy jet with similar initial conditions 
but different velocities. Figure [3] shows mass den¬ 
sity and velocity cuts in the r-z plane for jets with 
M = 10, 20, 30 and 150, respectively. In the four 
cases, the initial jet temperature is 10^ K and the 
ambient-to-jet density ratio is riambient/iijet = 0.1. 
The simulation was stopped at d ^ 8500 AU in 
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Fig. 3.— Two-dimensional mass density (upper 
half-panels) and temperature (lower half-panels) 
cuts in the r-z plane when the jet reaches 8000 AU 
from the beginning of the simulation. All the sim¬ 
ulations are for nambient/^jet = 0.1 but different 
Mach numbers. 

each case, for purposes of comparison. 

Last column in Table [T] shows the maximum 
temperature reached by the jet after impacting 
the dense cloud for each simulation. Temperatures 
similar to those observed in several X-ray emitting 
stellar jets (see Section [1]) are reached when the 
velocity of the jet is of the order of 700 kms“^ or 
higher. 

In Figure IH we summarize the results of our 
simulations. The top figure represents the tem¬ 
perature reached when the jet comes into the 
dense cloud {z ^ 120 AU in Figured]) as a func¬ 
tion of the ratio between the density of the cloud 
and the jet. The bottom panel shows the maxi¬ 
mum temperature as a function of the jet’s veloc¬ 
ity. No clear dependence of the temperature with 
^waii/'^jet is observed (top panel). At high den- 



Fig. 4.— Top: Maximum temperature of the jet 
after impacting the dense cloud for the different 
simulations, compared to the ratio between the jet 
and the dense cloud densities. Circle sizes are as in 
Figure O The dashed line indicates a temperature 
of 5 MK {kT = 0.43 keV). Filled circles are cases in 
which the jet reaches temperatures consistent with 
detected jets in X-rays. Bottom: Temperature 
distribution as a function of the velocity of the 
jet. 

sity ratios (uwaii/'^jet > 20), there is only a slight 



5 

















trend of increasing the jet’s temperature with in¬ 
creasing the density ratio. Filled circles in this 
figure are simulations with > 700 kms“^. On 
contrast, Figure IH bottom panel shows the tem¬ 
perature of the jet increases rapidly with increas¬ 
ing the jet’s velocity. Actually, the temperature 
reaches a plateau at ^ 1200 kms“^ where it re¬ 
mains approximately constant 10^ K) with in¬ 
creasing velocity. For > 700 kms“^, the jet 
reaches temperatures compatible with X-ray emis¬ 
sion (T > 10^ K). 

Summarizing, a high jet’s velocity is responsi¬ 
ble for the X-ray emission detected when the jet is 
lighter than the surrounding medium in the pro¬ 
posed scenario. 


4. HH 248: a candidate for the jet-cloud 
impact scenario 


The Herbig-Haro object HH 248 is situated to 
the South-East of the CTTS V615 Ori, at the 
base of the Horsehead Nebula (Barnar d 33) and to 


the South of the NGC 2023 nebula ( Malin et al. 


19871 ). The detection of an X-ray source at ^ 
15" from HH 248 and ^ 30" from V615 Ori 
(05^ 41”^ 25.7^, —0 2 ° 23^ 06.0") was reported by 
Lopez-Garcfa et al. ( 2013[ ) (their Src. 12). The 


authors did not find any optic al or IR counterpart 
down to 0.1 mJy. Src. 12 of iLopez-Garcfa et al ' 


(|2013[ ) is classified as an extended source by the 
XMM-Newton SAS source detection task, discard¬ 
ing an association with any possible highly ab¬ 
sorbed point-like source. The XMM-Newton ob¬ 
servation (ID 0112640201, revolution 419) was 
centered on the reflection nebula NGC 2023. The 
extended X-ray source was detected ^ 7' off-axis. 
Its proximity to the CTTS V615 Ori may be in¬ 
terpreted as they are related. However, an anal¬ 
ysis of the spectrum of this X-ray source gives 
A^h = 2 .3li;4 cm“^ (errors are at 90% 

confidence level) for the column density of the 
source and T r\j 10^ K as a lower limit. The 
A^h value is considerably larger tha n the one ob¬ 
tained bv lLopez-Garcfa et al.l ( 2013 ) for V615 Ori 
(A^h = 0.501o.o 4 cm“^). We conclude that 

V615 Ori is not the driving source of the extended 
X-ray emission detected. 


Assuming a distance of ^ 400 pc for the source. 
Lx = 3±2 X10^^ ergs“^ and EM « 2 x 10^^ cm“^. 
This analysis was performed by ourselves by fit¬ 



Energy [keV] 

Fig. 5.— Background-subtracted XMM/EPIC 
spectrum of the extended source including the 
model fit. 


ting an absorbed hot plasma model (jSmith et al 


200l[ ) to the observed EPIC/PN s pectrum using 


the X SPEC spectral fitting package (jArnan 


2004 ') with Cash statistics (ICa.shlll 9'^ . We notice 


that XSPEC treats both the source and the back¬ 
ground spectra separately when Cash statistics is 
used. The poor statistics of the source does not 
permit to study the properties of the X-ray emis¬ 
sion in more detail (see Eigure[5j). 

In the optical band, at an angular resolution 
of 0.3", HH 248 shows a knotty structure with 
five clear blobs. Eigure [6] shows an HST/WEC3 
image in the H band (Obs. ID 11548, obs. 
date 2010-10-19, exp. time 2496.170 s, PI: S. 
Megeath). A very faint source (WEC3 J05412764- 
0223143) is detected 3" to the South-East of the 
rad io source [RR94] NGC 2023 VLAl reported 
by iRodriguez fc ReipurthI ( 19941 ). at the posi¬ 
tion (a, 5) = (05^ 41“ 27.64", -02° 23' 14.3"). 
WEC3 J05412764-0223143 is connected to HH 248 
by a filamentary structure. WEC3 J05412764- 
0223143 is very likely the driving source of 
HH 248. Its magnitude in the E160W filter 
of the HST/WEC3 determined by us (see bel¬ 
low) is ^ 23.5 mag, corresponding to a flux 
Fjy = 0.46 ± 0.01 /iJy. Compared to the fluxes 
determined for the knots (see Table |2| this value 
is between two and three orders of magnitude 
lower. WEC3 J05412764-0223143 must be highly 
absorbed, with Ay > 13 mag according to its 
observed magnitude in the E160W filter. To de¬ 
termine the fluxes of WEC3 J05412764-0223143 
and the knots of HH 248, we performed aper- 
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05 : 41 : 28.8 05 : 41 : 27.6 05 : 41 : 26.4 05 : 41 : 25.2 


right ascension (deg) 

Fig. 6.— High-contrast HST/WFC3 H band im¬ 
age centered at HH 248 (Obs. ID 11548, PI: S. 
Megeath). Different high-density knots are iden¬ 
tified with capital letters. The circles mark radio 
sources identified in the literature. Their radius 
corresponds to the Icr positional error. The posi¬ 
tion of the extended X-ray source identified by us 
as the interaction of jet with the molecular cloud is 
also marked. The continuous line is the precession 
curve. 


Table 2: WFC3 F160W photometry of the knots 
and their probable driving source. 


Source 

AB 

(mag) 

F. 

(mJy) 

J05412764-0223143 

23.49 ± 0.18 

0.46 ± 0.01 

A 

18.21 ± 0.01 

58.96 ± 0.09 

B 

17.99 ± 0.01 

72.17 ± 0.10 

C 

18.95 ± 0.02 

29.88 ± 0.06 

D 

18.42 ± 0.01 

48.71 ± 0.08 

E 

16.56 ± 0.01 

270.16 ± 0.19 


Table 3: Initial physical conditions in the jet, the 
ambient (ISM) and the wall 


Parameter 

Value 

Jet 

Density 

Temperature 

Velocity 

500 cm“^ 

10^ K 

1.35 X 10^ kms“^ 

ISM 

Dentity 

Temperature®^ 

50 cm ^ 

10^ K 

Wall 

Dentity 

Temperature®^ 

5 X 10^ cm“^ 

10^ K 


ture photometry in the pipeline processed data 
were obtained from the Mikulski Archive for Space 
Telescopes (M AST) using the DAOPHOT package 
(Stetson 1987h . An aperture radius of 1" was used 
in each case to prevent contamination from other 
close knots, whose typical separation is ^ 2". We 
notice that a 1 "-radius extraction region encircles 
more than the 95% of each bright knot even for the 
case of knot E. Sky background value was fixed for 
subtraction. This value was estimated from the 
surroundings of the jet. An aperture correction 
was applied following the guidelines of the WFC3 
instrument 

The chain of knots in HH 248 follows a curve 
path typical of processed jets (e.g. 


Frank et al. 


20141 ). In Figur e [6] we show a simple fit of a preces¬ 


sion curve (e.g. Masque et al .1120121 ) to the knots 


^http: //www.stsci.edu/hst/wfc3/phot_zp_lbn. 


“Temperature is fixed to maintain pressure equilibrium 


and the IR source. The curve passes through the 
X-ray source. Hence, in the scenario proposed by 
us, the X-ray source is the bow shock produced by 
the encounter with a dense cloud of the jet flowing 
from WFC3 J05412764-0223143. 


Simulations for HH 248 


103 .^-3 


For our problem, we adopted a density of 5 x 
cm“^ for the wall as a lower limit, based of 
measurements of the density in the surroundings 
of the Horsehead Ne bula, where HH 248 is located 
( Habart et al. 2QQ5I ). Then we scaled the density 


of the jet to have a light jet traveling through the 
ISM (njet/'^iSM = 10) and a heavy jet through the 
wall (ujet/^waii = 0.1). Pressure equilibrium is 
imposed. As a consequence, the temperatures of 
the three regions are fixed. We assume an initial 
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Fig. 7.— Evolution of the jet’s intrinsic luminosity 
with time in the [0.3-10] keV energy band. The jet 
impacts onto the denser cloud at t = 1.1 x 10^ s. 


temperature for the jet of 10^ K. Thus, the tem¬ 
peratures of the ISM and the wall are, respectively, 
10^ and 10^ K. Table [3] summarizes the parame¬ 
ters of our model. Initially, the jet is defined as a 
homogeneous cylinder with size r ^ 100 AU and 
z ^ 1000 AU. The distance between the jet’s base 
and the dense cloud is kept as in Section O The ac¬ 
tual projected distance between the X-ray source 
and WFC3 J05412764-0223143 is 9500-12000 AU, 


assuming a distance from the Earth between 350 
and 450 pc as deduced from th e measurements for 
a Orionis (e.g. ICaballeroll2008[ ). 

The evolution of the jet’s luminosity with time 
is shown in Figure [71 X-ray luminos i ties we re de¬ 
rived as discussed in iBonito et al.l (2010b). In 
our space domain, the jet reaches the wall after 
^ 35 years. Figure [7| shows that the intrinsic lumi¬ 
nosity of the jet reaches ^ 10^^ ergs“^ after its im¬ 
pact against the molecular cloud. However, during 
the first stage of the jet’s evolution, while the jet 
crosses the (less dense) ISM, its X-ray luminosity 
remains at a lower value ( r\j 1029 gj.gg X- 

ray luminosity > 10^^ ergs“^ for the jet inside the 
high density molecular cloud is in agreement with 
the unabsorbed X-ray luminosity determined from 
the fit to the XMM-Newton datsU 



r[100 AU] 


Fig. 9.— Enlargement of the temperature (left 
panel) and density (right panel) maps derived 
from our simulations describing the jet/wall in¬ 
teraction (at 45 yr from the beginning of the sim¬ 
ulation). The complete movie of the evolution is 
also available as online material. 


Figure [8] shows the X-ray evolution map of the 
jet as it travels through the high-density molecular 
cloud. Each panel corresponds to a time delay 
of five years. The first panel is for 30 yr from 
the beginning of the model run. Figure [9] shows 
the temperature (left panel) and the density (right 
panel) maps describing the jet/wall interaction. 

The emission measure {EM) distribution near 
the end of the simulation (penultimate panel in 
Figure IHj) shows a small bump at T ^ 10^ K (see 
Figure Ho]). The mean value of the EM obtained 
from our simulations is compatible with the EM 
obtained from the fit to the X-ray data (^ 2 x 
10^^ cm“^). Note that the EM{T) at temperature 
< I MK shows a plateau analogous to that found 
by BOP07 in their simulations. This is mainly due 
to the radiative losses that are very efficient at T ^ 
10^ K and that cause a bump of EM{T) at T ^ 
10^ K. We have verified this point by comparing 
simulations either with or without the radiative 
losses. 


^The agreement with the observations is satisfactory be¬ 
cause we did not perform any fine-tuning of the initial 
conditions which is needed to accurately reproduce the ob¬ 
served luminosity. For instance, a better fit might be ob¬ 
tained by decreasing the ambient-to-jet density contrast 
(see BOP07). 


6. Discussion and conclusions 

In this work, we studied the scenario of a stellar 
jet impacting a molecular cloud from a HD mod¬ 
eling approach. In particular, we simulated su- 





















lE-07 2E-07 3E-07 4E-07 5E-07 6E-07 7E-07 8E-07 9E-07 

Fig. 8.— X-ray map of the jet after it reaches the high-density molecular cloud (wall). The period range is 
40-65 yr with a time binning of 5 yr. The pixel size if 1 AU. The color scale is cnt s“^ pixel“^. 


personic stellar jets moving through ISM medium 
that is less dense than the jet encountering a dense 
molecular cloud. Hence, the jet becomes under- 
dense inside the molecular cloud. From previ¬ 
ous results (e.g. BOP07), X-ray emission is ex¬ 
pected when the jet is lighter than the surround¬ 
ing medium, while it is not when the jet is heavier 
than the ambient. 

Our simulations explore the parameter space: 
jet’s velocity, ISM density and molecular cloud 
density. The scenario of a light jet that impacts a 
dense cloud may be typical of some stellar jets ob¬ 
served in star-forming regions. In particular, this 
scenario is possibly occurring for HH 248. Our 
main goal is to investigate whether such an im¬ 
pact may produce X-ray emission from the jet. 

Our results from the simulations indicate that 
jets with v^et > 700 kms“^ can reach temperatures 
of several MK after impacting the dense molecular 
cloud. This result is in agreement with detection 
of X-ray emission from the jet. In contrast, low 
velocity jets do not reach such temperatures. We 
find a clear trend of increasing temperature after 
impact with increasing the jet’s velocity, indepen¬ 


dently of the density ratio between the jet and the 
dense molecular cloud. 


We used these results to investigate the possible 
origin and evolution of HH 248, including the pos¬ 
sibility of X-ray emission after encountering the 
high density region at the base of the Horsehead 
Nebula. From our multiwavelength analysis of the 
HH object, we reveal the probable driving source 
of the jet, which is a highly absorbed point-like 
source detected only in a deep H-band image with 
the HST/WFC3. From the WFC3 image, we iden¬ 
tify five dense knots that are not fully aligned with 
WFC3 J05412764-0223143. We used them to de¬ 


termine the parameters of the precession of the 
jet associated with HH 248. The fitted preces¬ 
sion curve pass through the ext e nded X-ray source 


Sion curve pass tmougn tne ext e nded A-ray source 
named bv Ihopez- Garcia et al. I (l2ni.^ as Src. 12. 


We conclude that this X-ray source is the front 
shock of the jet that created HH 248. 


From our X-ray data analysis, we discard the 
possibility that the extended X-ray source were 
part of the PSF of the classical T Tauri star 
V615 Ori, also detected in the XMM-Newton as 
^ 30". In particular, the column density de- 
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Fig. 10.— Emission measure distribution {EM) 
as a function of temperature (T) for the penulti¬ 
mate panel of Figure [HI at a time of 25 yr after 
the interaction of the jet with the molecular cloud 
(60 yr from the beginning of our simulation). 


termined for this source is one order of magni- 
tude higher than that dete rmined for V615 Ori by 


Lopez-Garcfa et ah! (|2013[ ) . We propose a scenario 


for HH 248 in which the jet, which moves through 
the ISM, reaches a region of much higher density. 
This region may be identified with the ionization 
front IC 434 that runs south from the NGC 2024 
star-forming region and contains the Horsehead 
Nebula (Barnard 33). In this scenario, the jet, 
which is heavier than the ISM but lighter than the 
molecular cloud, emits in X-rays when it reaches 
the denser cloud. Hence, the X-ray emission from 
the jet would be produced far from the jet base, 
contrarily to other c ases (see Sect i on [T]). In the 
classification given bv iBonito et al.l (2010b), lumi¬ 
nous X-ray jets detected at large distances from 
their driving source seem to be associated with 
high-mass stars. In our case, the location of the 
X-ray emission from the jet would not be related 
to the mass of the driving source but to the change 
in the conditions of the ISM. This is an excellent 
test case to check the idea that detectable X-ray 
emission from stellar jets is produced if the jet 
travels through a denser environment (BOP07). 

From our results, we conclude that X-ray emis¬ 
sion from the impact of a stellar jet with a dense 
cloud is feasible. Deep X-ray observations of HH 
objects in nearby star forming-regions may reveal 
other cases similar to HH 248. New X-ray ob¬ 
servations are needed to determine precisely some 


parameters of the X-ray counterpart of the jet. 
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at the University of Chicago. HST data presented 
in this paper were obtained from the Mikulski 
Archive for Space Telescopes (MAST). STScI is 
operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA con¬ 
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